The pattern of charge ordering in quasi-one-dimensional organic charge transfer solids 
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We examine two recently proposed models of charge ordering (CO) in the nominally ^-filled, quasi- 
one-dimensional (ID) organic charge transfer solids (CTS). The two models are characterized by site 
charge density "cartoons" ...1010... and ...1100..., respectively. We use the Peierls-extended Hubbard 
model to incorporate both electron-electron (e-e) and electron-phonon (e-ph) interactions. We first 
compare the results, for the purely electronic Hamiltonian, of exact many-body calculations with 
those of Hartree-Fock (HF) mean field theory. We find that HF gives qualitatively and quantitatively 
incorrect values for the critical nearest-neighbor Coulomb repulsion (V c ) necessary for ...1010... order 
to become the ground state. Second, we establish that spin-Peierls (SP) order can occur in either 
the ...1100... and ...1010... states and calculate the phase diagram including both on-site and 
intra-site e-ph interactions. Third, we discuss the expected temperature dependence of the CO 
and metal-insulator (MI) transitions for both ...1010... and ...1100... CO states. Finally, we show 
that experimental observations clearly indicate the ...1100... CO in the 1:2 anionic CTS and the 
(TMTSF) 2 X materials, while the results for (TMTTF) 2 X with narrower one-electron bandwidths 
are more ambiguous, likely because the nearest neighbor Coulomb interaction in these materials is 
near V c . 

PACS numbers: 71.30.+h, 71.45.Lr, 75.30.Fv, 74.70.Kn 
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I. INTRODUCTION 



Recent experimentsU'BM'BlMlid'BBme 
showing clear evidence for the existence of charge-order 
(CO) in several 2:1 cationic organic charge transfer 
solids CCTSQ iMr^H i t 1J iWr a fr 5 ^ considerable theoretical 
interestBBM»BBE3BM As these materials are 
nominally i-filled (one electron or hole per two sites) and 
involve electron-electron (e-e) interactions (including 
nearest- neighbor Coulomb repulsion V), the apparently 
obvious charge ordering in the quasi-lD CTS systems is 
the Wignep-jpystal-like "...1010..."|-*tate--a»d^umerous 
theoreticaEKa and experimentaBllSllS'Htata studies 
have argued for this possibility. The site occupancies 
T' and '0' actually correspond to 0.5 + e and 0.5 - e, 
respectively, with e an important measurable quantity. 
The ...1010... CO corresponds to a 4k F CDW (spatial 
period 2), where kp = 7r/2a is the Fermi wavevector in 
the absence of Coulomb interactions, with a the riaitieft 
constant. On. the.pther hand, in previous work weOiiiatS 
and othersE3-E2lE3LJ have shown that for realistic values 
of V and when electron-phonon (e-ph) interactions are 
included, a "...1100..." charge ordering can become 
the ground state, and, importantly, that this CO state 
can explain the "mysterious" states observed in several 
CTSlirau, in which coexisting charge and spin density 
waves of the same periodicity occur. The ...1100... CO 
corresponds to a 2kp CDW (spatial period 4) and has 
a co-operative coexistence with the 2k F BOW or a 
mixed 2kp + 4kp BOW, which is why pwc-tpjined it a 
"Bond-Charge Density Wave" (BCDW)BBE2. In the 
present article, we examine critically the theoretical and 
experimental evidence for each of these two possible 



charge orderings in the quasi-lD organic CTS. Before 
proceeding to our analysis, it is useful to formulate in 
a precise manner several key questions related to the 
charge ordering in the organic CTS. We focus here on 
four such questions. 

First, the ...1010... CO ground state requires both 
strong on-site Coulomb repulsion (U » to within 
the extended Hubbard model) and nearest-neighbor 
Coulomb repulsion (V within the extended Hubbard 
model) to be greater than a critical value, V C (U). For 
any comparison to real materials determining the correct 
value V C (U) is essential. Both Hartree-Fock (HF) and 
exact many-body methods have been used to estimate 
V C (U). How reliable are the HF estimates quantitatively, 
and do they exhibit the correct qualitative trends as U 
varies? 

Second, charge ordering is only one of the phenom- 
ena observed in organic CTS. At lower temperatures, 
broken symmetry states involving spin ordering — spin- 
Peierls (SP) and spin density wave (SDW) states — are 
typically observed, to say nothing of superconductivity. 
Many of the quasi- ID CTS materials that show CO also 
feature S,P, .ground states at low temperature. Our previ- 



ous workllZrEatJ establishing the BCDW state proves that 
the 1100 CO is consistent with SP. To our knowledge, 
none of the t he o re tical studies proposing the existence of 
...1010... COBEl has investigated whether a SP phase 
can occur within the ...1010.... CDW at lower temper, 
atures. Indeed, an early quantum Monte Carlo studyEZl 
can be taken to suggest that a SP phase could not coex- 
ist with the ...1010... CDW. Is 1010 CO consistent with 
a SP phase at low temperatures, and if so, under what 
conditions? 

Third, any successful theory of charge ordering within 
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FIG. 1: Phase diagram of the ID extended Hubbard model at 
quarter-filling for positive values of U and V. The 4k^ CDW 
exists only in the upper right hand corner, above the value 
of V C (U) determined by the two strong coupling curves (solid 
line: U — > oo; dashed line: V — * oo) and the quantum Monte 
Carlo data points. Dashed-dotted line: Hartree-Fock data. 



the organic CTS must explain not only the mechanism 
and pattern for CO, but also the temperature scale at 
which it occurs. In particular, any complete theory of CO 
in these materials must be able to explain the existence of 
different temperatures for the metal-insulator (MI), CO, 
and SP transitions. Does either of the theoretical models 
for CO explain the observed sequence of transitions? 

Fourth, the ultimate arbiter in the ...1010... vs. 
...1100... debate is experiment, and in this regard it is 
clearly desirable to have as broad a range of probes and 
materials as possible. Both the 2:1 cationic CTS (with 
holes as carriers) and the 1:2 anionic CTS (with electrons 
as carriers) have been studied extensively using similar 
theoretical and experimental approaches. What is the 
experimental evidence for ...1010... and for ...1100... in 
both these classes of materials? 

We address these questions in the remainder of this 
paper. In Section || we introduce the the electronic ex- 
tended Hubbard model and compare exact to mean field 
calculations. Here we establish that HF theory exagger- 
ates the role of the ...1010... CO and that the realistic 
range of parameter values for which this CO is obtained is 
rather narrow. In Section HI we introduce e-ph couplings 
and show that it is possible for a spin Peierls transition 
to occur with the 



...1010... CO phase In Section |I_V|, 
we discuss the theoretical expectations for the sequence 
of temperatures at which the various transitions-metal- 
insulator (MI), charge order (CO) and spin-Peierls (SP)- 
occur, establishing that within the strictly ID ...1100... 
CO model, Tmi > Too = Tgp, whereas within the 



strictly ID ...1010... CO model, T MI = T co > T SP , 
In Section ^ we compare our theoretical results with ex- 
perimental data on a wide variety of CTS, both cationic 
and anionic. These comparisons establish that in both 
the anionic 1:2 CTS and the (TMTSF) 2 X the dominant 
CO pattern appears to be ...1100..., which in the latter 
case coexists with a SDW (because of two-dimensional 
effectsEI) rather than a SP phase. For the (TMTTF) 2 X 
CTS, the situation appears less certain, but ...1010... or- 
der appears to be favored. We propose specific experi- 
ments to distinguish definitively between the two differ- 
ent CO patterns. We also show in Section IV that the 
experimentally observed temperatures of the MI, CO, 
and SP transitions are not fully consistent with either 
...1010... or ...1100... models in their strictly ID forms. 
Finally, in Section we summarize our results and list 
some open problems for future research. 



II. THE ID EXTENDED HUBBARD 
HAMILTONIAN 

We first consider the ID extended Hubbard Hamilto- 
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In the above, j is a site index, rij — rij^ + and a 
is spin. We focus on the i-filled case, with the average 
number of electrons (or holes) per site, p = \ . In addition 
to the purely electronic terms in Eq. ([!]), we will also 
consider e-ph couplings explicitly in section P4 

That for U — > oo there exists a critical value (V c ) 
of V for the appearance of the 4k F CDW (...1010...) 



detailed studics-pf V c for finite U have appeared in re- 
Nonetheless, it appears that the im- 



within the ^-filled ID "pure" extended Hubbard model 
in Eq. ([!]) has been known for decadesa, and several 

cent yearsrfff 
plications of these theoretical results have been insuffi- 
ciently appreciated in the recent literature on organic 
CTS. Thus, we believe these results and the underly- 
ing arguments bear revisiting briefly here. In the limit 
U —> oo, the ^-filled band of spin-i electrons becomes 
equivalent to a half-filled band of spinless fermionsLStS, 
which in the presence of V can be mapped using a 
Jordan- Wigner transformation to the Heisenberg XXZ 
chain. This rigorous mapping and the exact solution for 
the XXZ model establish that for V > 2to, the system is 
in the (gapped) "Ising-Heisenberg" phase, which in terms 
of the original electronic model corresponds to the 4kp 
CDW order ...1010... For V < 2to, the spin system is in 
the "Heisenberg-XY" phase and there is no gap, which 
in the pure extended Hubbard model corresponds to a 
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FIG. 2: (a) Charge disproportiation An vs. V within the 
HF approximation for U = 6io. Open (filled) symbols are for 
iV = 16 (iV=32) rings, (b) V c within the HF approximation. 
Lines are guides to the eye. 



Luttinger liquid (LL) with no charge order. Hence for 
U — > oo, the critical value of V is V c (oo) = 2to. 

For finite U, the value of V c at finite U has been 
calculated i) within strong coupling perturbation the- 
ory around U — > oo, where a second-order calculation 
yields^ V C {U) ~ 2t + + ii) within strong- 

coupling theory in the limit V — > .tra, where a sim- 
ilar second-order calculation showscSEl that for finite 

V the 4kp CDW phase boundary occurs at U C (V) — 
4t + yj- + • • • ; and iii) for intermediate values of U and 

V by both exact.|diagonalizatiorO and quantum Monte 
Carlo methodsEiM In Fig. | we show the resulting 
boundary of the 4kp CDW phase in the L 7 , V plane; in 
this figure, the solid and dashed curves represent the two 
strong coupling expansions and the points with error bars 
represent our new determination of the boundary at in- 
termediate coupling using a quantum Monte Carlo ap- 
proach described in Ref. |3l|. 

The relevance of presenting the true phase boundary 
in the U, V plane (as determined by accurate many- 
body methods) becomes clear when one considers the 
results obtained in a number of recent studies of the 
i-filled ID extended Hubbard model using mean-field 
approximationsEJc3c3. Using the HF approximation for 
the U interaction, and Hartree for the V interaction, as is 
done in references E20L leads to drastically reduced values 
of V c compared to the true results: for U = 5to, the mean- 
field V C MF = OAtjB, as compared to V c = (3.15±0.07)f 
expected from Fig. |l|. In Fig. ^, we present results for 
our own unrestricted HF calculations on finite size lat- 
tices. In these calculations, both U and V interactions 
were treated within the full HF approximation, and no 
assumptions are made about the periodicity of the solu- 



tion. The value of V c is readily determined by the point 
at which An, the difference between charge densities on 
adjacent sites, becomes nonzero (see Fig. 0(a)). As seen 
in Fig. 0(a), finite size effects are quite small, and V c 
for the 64 site lattice did not differ significantly from 32 
and 16 sites. From Fig. ||(b), we see that contrary to the 
known exact results for the ID extended Hubbard model, 
within the HF approximation (i) the ...1010... CDW does 
occur for U < ito, and (ii) the HF V c increases with U. 
For direct comparison of exact and HF results, the data 
of Fig. |2|(b) are plotted as the dashed-dotted line in Fig. |l| 
and clearly show that the mean-field approximations in- 
correctly increase the parameter regime in which ...1010... 
CO is expected in the extended Hubbard model. 

III. CHARGE AND SPIN ORDER INCLUDING 
PHONONS IN ID 

To answer the second question posed in the introduc- 
tion, we explore systematically the quantitative effects of 
adding e-ph interactions to the extended Hubbard model 
in Eq. (|l|), resulting in a "Peierls-extended Hubbard" 
(PEH) model. Our results establish that a SP transi- 
tion can occur within the ...1010... CO state and clarify 
the regions of e-ph parameter space in which the two 
competing forms of charge order exist. Our PEH model 
follows from replacing the one-electron part Hq in Eq. (Q) 
by 

Hq = - Y^[t0 ~ "A; />',,. | +P^2 VjUj 

In the above Aj = Uj+i — Uj, where Uj is the dis- 
placement of the jth atom from equilibrium. The am- 
plitude of the internal molecular vibration is given by 
Vj. a and (3 are the intersite and intrasite e-ph coupling 
constants, respectively, and K\ and are the corre- 
sponding spring constants. The kinetic energy operator 
Bj,j+i,a = c j+i, < r c j> + c i<r c 3+ 1 ' <r ' We determine (numeri- 
cally) self-consistent ground state solutions for the above 
PEH model as functions of the parameters, measuring 
charge density (rij) and bond order Bj,j+i,o-)- 

We ,-s.tar.t M^jith a summary of previous relevant 
resultsOEJ'EiJoEj concerning the SP transition in the 
i-filled band. These prior investigations have estab- 
lished that for V < V c and for non-zero a or f3 (or 
both), the dominant broken symmetry below the metal- 
insulator but above the insulator-insulator SP transition 
is the dimerized 4ki? BOW with uniform site charges 
(Fig. §(a)). The SP transition here is a second dimer- 
ization of the dimerized lattice, leading to the BCDW 
state, which can be considered as a superposition of the 
2k^ and 4kp BOW's, and which is accompanied by the 
...1100... CO (Fig. |(b)). The bond distortion pattern in 
the BCDW is SW'SW, where S is a strong bond, and W 
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FIG. 3: The competing density wave ground states in the 
Peierls-extended Hubbard Model at i-filling. Grey, black and 
white circles correspond to site charges of 0.5, 0.5 + e and 0.5 
- e, respectively, (a) The BOW state, with dimerized 
bond orders SWSW and uniform site charge densities, (b) 
The BCDW state, with bond orders SW'SW and the accom- 
panying CO ...1100... (c) The 4k F CDW with ...1010... CO 
and uniform bond order, (d) The 4fcf CDW-SP state, with 
...1010... charge order and bond orders SSWW. The charges 
on the "unoccupied" sites in (d) are slightly different (see 
Fig. |). 



and W are weaker bonds with W' > W. In contrast, the 
bond distortions within the 4kp CDW-SP phase is of the 
form SSWW. The competing CO is the ...1010... 4k F 
CDW, which has uniform bond orders at high temper- 
atures (Fig. |(c)). The SP phase within the 4k F CDW 
phase, whose existence has not previously been studied or 
established, would correspond to Fig. |3|(d). We demon- 
strate below the existence and nature of this phase, which 
we hereafter refer to as the 4k F CDW-SP phase. Within 
the 4kp CDW-SP phase there occur alternate strong and 
weak bonds between the "occupied" sites that are ac- 
tually second neighbors (see Fig. ||(d)). Previous work 
did not find this 4kp CDW-SP— nphaae either because of 
focusing primarily on V < V c iIi&^B, or because of in- 
cluding only one, but not simultaneously both of the e-ph 
couplings a and /Xj'Ej. In this section, we present new 
results for the i-filled band for both V <V C and V > V c , 
including both e-ph couplings a and 8. 

We have performed fully self-consistent exact diag- 
onalizations of finite-sized periodic rings. The self- 
consistency equations derived from 
are 
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We present results in terms of dimensionless e-ph cou- 
plings X a — a 2 /(Kit ) and A/3 = /9 2 /(i^ 2 io)- We have 
obtained results for systems of sizes N=8, 12, and 16. In 
Fig. |J we present a ground state phase diagram (iV=16) 
in the X a — A/3 plane, showing regions in which the various 
broken symmetry states dominate. In all cases, U = 8to; 
panels (a), (b) , and (c) correspond to V = 2t ,3t , and 
4io, respectively. 

Consider panel (a) with V = 2t . Since U — 8t , 
this value of V lies below V c for the infinite system (see 
Fig. 0). For this 16-site system, as we increase X a from 
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FIG. 4: The phase diagram in the X a — Xp plane for U = 8to 
and (a) V = 2t ; (b) V = 3t ; and (c) V = 4t for a 16 
site periodic ring. The boundaries between the the various 
phases, drawn as continuous lines in the figures, are actually 
determined numerically with an accuracy of ~ 0.05 in X a and 
Xp. In the regions labeled "unphysical" the weakest bonds 
have negative bond orders due the excessive lattice distortion 
that occurs for very large X a . 



zero (for fixed small Xp), there is initially no distortion. 
For A Q > 0.7, the ground state becomes the 4kp BOW 
shown in Fig. 0(a); while the site charge densities re- 
main uniform (n=0.5), the bond-order is now inhomo- 
geneous and has the form SWSW. For X a > 0.9, the 
BCDW discussed above and depicted in Fig. |||(b) be- 
comes the ground state. Hence the bond-order pattern 
now is SW'SW and the charge order has the ...1100... 
pattern. Next consider increasing Xp from zero at fixed 
small X a . Again the state is initially undistorted, but 
when A/3 reaches a critical value, the Akp CDW with 
...1010... charge order and uniform bond order becomes 
the ground state. Importantly, we see that a SP dis- 
tortion does occur within the ...1010... charge-ordered 
region, but the resulting Akp CDW-SP state (shown in 
Fig. ||(d)) occurs only for fairly large values of both X a 



and A/3 when U = 8to and V = 2t , at least for this 16- 
site system. Finally, we note that in the region labeled 
unphysical in Fig. || the value of A Q has become so large 
that self-consistency drives the weakest hopping integral 
negative, rendering the results unphysical in this context. 

Panels (b) and (c) of Fig. |] show several important 
trends as V is increased through and beyond the value 
of V c (8to) for an infinite system. First, the sizes of the 
undistorted, 4fcp BOW and BCDW regions decrease con- 
siderably with increasing V , and move to larger values of 
A a , showing that it takes stronger e-ph coupling to over- 
come the "natural" tendency toward ...1010... CO at 
large V. For the same reason, the sizes of the Akp CDW 
and Akp CDW-SP regions increase considerably with V 
and occur for smaller values of \p. Second, the existence 
of the CDW-SP ground state is very robust — V does 
not appear to suppress this SP order. Third, Fig. || il- 
lustrates one clear qualitative difference between the two 
SP phases in this finite size system. Provided that X a is 
greater than a critical value depending on V (A Qc (V)), 
the BCDW occurs for all \p > 0. In contrast, occur- 
rence of the Akp CDW-SP state requires both A Q and 
\p to be non-zero. Note however that as V is increased, 
the minimum value of A/3 (call it Xp c (V)), for which the 
Akp CDW-SP state exists decreases rapidly towards zero, 
while the value of A Q necessary to produce this phase at 
A/3 c (V) increases slightly. 

These intriguing finite size results beg the question of 
what happens to the system in the thermodynamic limit, 
N — > oo. Unfortunately, our limited data do not allow us 
to perform a reliable finite size scaling analysis, and we 
can thus provide only tentative, partial answers at this 
stage. We believe that in the N — > oo limit the undis- 
torted and Akp BOW phases disappear; our evidence for 
this comes from our earlier studies of the BCDW (hence 
for V < V c ) on long open chainsEj, which showed that the 
center of these chains distorts naturally into a BCDW but 
which showed no sign of a pure Akp BOW. Similar stud- 
ies for V > V c have thus far proven inconclusive, and the 
N — > oo behavior in this regime remains an important 
open question. 

Our self-consistent exact diagonalizations give results 
for site charges, bond orders, and hopping integrals for all 
the phases, homogeneous and inhomogeneous. In Fig. ^ 
we present the actual site charge densities and bond or- 
ders for U = 8to, V = Mq, and A/3=0.6. This plot corre- 
sponds to a vertical "slice" through Fig. |](b), and shows 
the 4fc f -CDW, 4fc F -CDW-SP, and BCDW phases. Sev- 
eral points are suggested by the data in Fig. [| 

First, in the Akp CDW-SP state, the sites with smaller 
charge density are no longer equivalent — there is now (as 
expected from the bond distortion) a small charge differ- 
ence between these two sites. Thus, the Akp CDW-SP 
state is characterized by two different bond orders and 
three different site charges, as opposed to the two differ- 
ent charges and three different bond orders of the BCDW 
state, which are also seen clearly in Fig. |[ We shall re- 
turn to this point momentarily. Second, the quantitative 
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FIG. 5: Different possible values for charge densities and bond 
orders for U = 8t , V = 3t , and X = 0.6 (see Fig. gb)). 
Circles (diamonds) are data for N—16 (N—12) site rings. 



values of the site charge differences, Art, in the BCDW 
state are quite small, whereas the An in the 4/cf CDW 
and 4/cf CDW-SP states are considerably larger. We 
have not been able to find, in our finite system calcula- 
tions, any BCDW state with large An, nor any AUf CDW 
or 4kp CDW-SP states will small An. Third, and simi- 
larly, the bond order differences in the BCDW are always 
quite large, whereas the bond order differences in the AUf 
CDW-SP phase are always quite small. Again, these in- 
triguing results for finite size systems beg the question 
of the behavior for N — > oo, and again we cannot pro- 
vide definitive statements because our data are not suf- 
ficient to allow a reliable finite size scaling analysis. One 
argument suggesting that the large An may not be an 
inevitable consequence of the ...fOfO... CO comes from 
the U — * oo limit and the mapping to the Heisenberg- 
Ising model; translating the exact results into fermion 
language shows that charge difference in the ...1010... 
CO phase is given by An ~ exp — (y^ 2 t ) > so that it 
starts from zero-and remains small for a large region of 
V > V c — 2t c3- Clearly, further study of these po- 
tentially significant differences in the N — > oo limit is 
required. 

Beyond the quantitative distinctions, however, there is 
the clear qualitative distinction between the BCDW and 
the Akp CDW-SP states: namely, the former always has 
two different site charges and three different bond orders, 
whereas the latter has three different site charges and two 
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FIG. 6: Schematic showing expected temperature dependence 
of (a) ...1010... and (b) ...1100... CO models. T S p and Tmi 
indicate SP and MI transition temperatures respectively. The 
shaded region indicates temperature range where CO is found. 



In the ...1100... CO— model, one must consider a 
quadruple-well potential!^. At the highest tempera- 
tures there are four equivalent configurations, ...1100..., 
...0011..., ...0110..., and ...1001..., again giving a metallic 
state when all four have equal weights. The high tem- 
perature transition corresponds to a symmetry breaking 
that leaves ...1100... and ...0011... with equal weight, 
resulting in no net CO, but an insulating state. Then 
at lower temperature, the ...1100.../... 0011... symmetry 
breaks, resulting in CO. As the ...1100... CO is a co- 
operative BCDW state, CO and bond distortion occurs 
simultaneously^. Thus within the ...1100... model (sec 
Fig. ^|(b)), one again expects two transitions, a high tem- 
perature MI transition, followed by a low temperature 
combined CO/SP transition. 



different bond orders. This qualitative distinction will 
persist in the large system limit and is in principle ac- 
cessible to experimental techniques sensitive to the local 
charge or bond order environment. For instance, if NMR 
measurements can be made below the spin Peierls tran- 
sition temperature (T$p) in a given CTS, they should be 
able to distinguish cleanly between the BCDW and the 
4fcp-CDW-SP states. We will return to this and related 
points in our comparisons to experiment. 



IV. 



TEMPERATURE DEPENDENCE OF CO 
TRANSITIONS 



Our calculations in Sections III and ||, as well as 
virtually all other theoretical studies, have been lim- 
ited to ground state results only. However, the issue of 
the temperature dependence of the CO and MI transi- 
tions is quite important in applying the results of ID 
model calculations to experimental systems. In this Sec- 
tion, we discuss the temperature dependence expected 
for ...1010... and ...1100... models of CO. 

As the ...1010... CO is driven by V, it is expected to 
give CO already at quite high temperatures. The series 
of transitions can be understood by recognizing that the 
...1010... CO has two equivalent configurations, ...1010... 
and ...0101..., corresponding to a double- well potential. 
For temperatures above the CO transition temperature, 
both possible configurations ...1010... and ...0101... have 
equal weight in the partition function, resulting in no net 
CO and a metallic state. The CO transition then corre- 
sponds to a symmetry breaking between ...1010... and 
...0101..., with the result that CO and MI transitions (at 
Tqo an d Tmi respectively) happen at the same temper- 
ature. As we have shown in Section III, a SP distortion 
can occur within the ...1010... CO phase, but its scale 
is set by the e-ph couplings and is expected to be well 
below the electronically driven Tmi = Tco ■ Thus within 
the ...1010... model (see Fig. ||(a)), one expects two tran- 
sitions, a combined high temperature CO/MI transition 
followed by a low temperature SP transition. 



V. INTERPRETATION OF EXPERIMENTS 

If detailed comparisons of models to real systems are 
to be attempted, accurate estimates of the e-e and e- 
ph parameters are required. Determination of U and 

V from quantum chemistry calculations is difficult, and 
in general such calculations give values that are too 
largeEJ. Fortunately, there is by now general agree- 
ment on the magnitudes of to in the different materi- 
als, and reasonable values of U and V can be estimated 
from the experiments. We consider primarily materials 
based on the TCNQ, TMTSF, and TMTSF molecules. 
It is generally accepted that the values of to in the 
cationic TMTTF/TMTSF and anionic TCNQ CTS lie 
between 0.1 - 0.25 eV. We then agree with previous in- 
vestigators that U/t in (TMTTF) 2 X raiLges from 7 - 
12E3E9 and is less than 6 in (TM^F^e^^ESEZl. Based on 
comparing relative 4kp responsesE3'E2le3 (a signature of 
large U/to), one may conclude that for TCNQ systems 
U/to is intermediate between the TMTTF and TMTSF 
limits. There are far fewer estimates of V/to' values 
of 2.8 for (TMTTF) 2 X_.and 2.0 for (TMTSF) 2 X have 
been proposed by MilaE£l. We note that whether or not 

V > V C (U) can be determined from the experimental pat- 
tern of the CO, if it is known. On the other hand, given 
the relatively small U/to in (TMTSF) 2 X, it is unlikely 
that V > V C (U) in this system (recalling the additional 
condition that V < U/2, based on general considera- 
tions of the nature of realistic Coulomb interactions.). 
To summarize, in the quasi- ID materials, we believe 



that (U/t ) 



TMTTF 



> (U/h) 



TCNQ 



> (U/to) 



TMTSF 



For the e-ph couplings, the dimensionless coupling con- 
stants allK.uta and j3/Kito are even more difficult to 
estimateoocj. For our overall purposes, it suffices to 
note that there is evidence for both nonzero a and f3. Fi- 
nally, although in this paper we do not explicitly model 
the effects of transverse coupling of the ID chains, (see 
reference |l| for a quantitative discussion of these effects) , 
we shall comment on these higher dimensional effects in 
some of the experimental discussion below. In this regard 
it is important to recall that the electronic anisotropy and 
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effective e-e interaction strengths can vary independently 
among the CTS. The relative strength of the e-e interac- 
tions is measured by the ratios U/to and V/to, whereas 
the effective dimensionality is measured by t±/to). These 
three ratios vary considerably among the materials. In 
particular, the TMTTF materials are both highly ID and 
very strongly correlated, whereas the TMTSF materials 
are less ID and also less strongly correlated. 

Let us turn now to discuss the experimental situations 
with individual classes of materials. 



A. 1:2 anionic CTS 

As mentioned in the Introduction, we believe that a 
complete understanding of the CO phenomenon in the 7- 
filled band CTS requires that both cationic and anionic 
CTS be examined on the same footing. The availability 
of large single crystals in the case of the anionic sys- 
tems makes these particularly attractive. On the other 
hand, the anionic CTS are much older and many dif- 
ferent contradictory claims exist in the literature. We 
therefore focus only on the materials for which clear ev- 
idence for CO has been found. The systems we consider 
are MEM(TCNQ) 2 , TEA(TCNQ) 2 , (DMe-DCNQI) 2 Ag, 
and (DI-DCNQI)2Ag. Preliminary discussions of the CO 
in MEM(TCNQ) 2 and TEA(TCNQ) 2 were also given in 
reference pj[ 

The most direct evidence for the ...1100... CO has 
been found in MEM(TCNQ) 2 and TEA(TCNQ) 2 . The 
SP transition in MEM(TCNQ) 2 occurs at 17.4 K, and 
a very careful neutron diffraction- study was performed 
for deuterated samples at 6 KcJ. This work clearly 
established the SW'SW bond distortion pattern asso- 
ciated with the ...1100... BCDW for T < T SP in 
MEM(TCNQ) 2 . Similarly, the charge and bond modu- 
lation pajtt|pnJtjit,TEA(TCNQi^,have been investigated 
by X-rayrirTrlr 8 ! and neutronEJ measurements. Below 
the structural transition at 210K, the amplitude of the 
SW'SW distortion increases rapidly, and ...1100... CO 
was foundocZI. 

The pattern of the CO in (DI-DCNQI) 2 Ag (but not 
in (DMe-DCNQI) 2 Ag), was until recently controversial. 
The (DI-DCNQI) 2 Ag system was originally claimed to 
show the ...1010... CO from NMR experimentsO. It was 
pointed out by the present authors that some of the ex- 
perimental results, especially the appearance of the CO 
at a temperature lower than the metal-insulator tran- 
sition-, (see Section |iv| ), perhaps indicated the ...1100... 
ordem. Initial optical measurements had detected evi- 
dence for a 4k F BOW without CO already at room tem- 
perature, which had suggested that a ...1100... CO ac- 
companies the second lattice distortion that occurs at 
lower temperatureEl. A first X-ray structural study found 
moderately strong satellite intensities in the 4k F phase 
(10~ 3 pi the Bragg reflection), accompanied by high I 
indicesE^.j-This was interpreted to indicate a displacive 
transitioned, suggesting (though not necessarily proving) 



that the 4k F phase is a BOW and not a CDW. This 
would again suggest that the CO that occurs at low tem- 
perature is the ...1100... CO. A later X-ray experiment 
included both (DI-DCNQI) 2 Ag apd the structurally re- 
lated material (DMe-DCNQI) 2 AgEil. The strong 4k F re- 
flection in the latter material (10 -3 - 10~ 2 of the Bragg 
reflection), along with the occurrence of 2k F reflections, 
indicates that the 4k F phase is a BOW (as agreed upon 
also by the authors of Ref. |5l]) . This would again indicate 
the ...1100... CO in the 2k F phase. No 2k F reflections 
were found in (DI-DCNQI) 2 Ag, and the index depen- 
dence of the 4k F reflections were also different from that 
in (DMe-DCNQI) 2 AgLli Based on this, the authors of 
Ref. |l| suggested that unlike in (DMe-DCNQI) 2 Ag, the 
4k F order in (DI-DCNQI) 2 Ag is a CDW, although their 
calculated satellite intensities, with the assumption of a 
...1010... CO accompanied by molecular distortions, were 
still lower than that of the observed intensities. The oc- 
currence of a 4k F BOW in (DMe-DCNQI) 2 Ag and a 4k F 
CDW in (DI-DCNQI) 2 Ag, though not impossible, would 
have been rather mysterious. This mystery has very te.- 
cently been resolved by Menghetti and collaboratora£3. 
The authors have analyzed the temperature-dependent 
vibronic and vibrational infrared absorptions in the (DI- 
DCNQI) 2 Ag system in great detail and have concluded 
that the low temperature phase is a 2k F bond tetramer- 
ized phase with the ...1100... CO pattern. Indeed, the 
authors also find that A n in this system is rather small, 
in agreement with the predictions of our finite size cal- 
culations in Section [II. 



To conclude, the CO pattern in all 1:2 anionic CTS 
for which detailed experimental data are available is 
...1100.... 



B. (TMTTF) 2 Br and (TMTSF) 2 X 

(TMTTF) 2 Br exhibits behavior similar to that of 
(TMTSF) 2 X, in that the. low temperature phase here is 
the mixed CDW-SDWHtl and different from the other 
TMTTF-based materials. We therefore discuss this ma- 
terial together with the (TMTSF) 2 X. There is strong ev- 
idence for the ...1100... CO in this class of materials, as 
we now discuss. 

First, we have pointed out in the above that U/to and 
V/to in TMTSF are very likely smaller than that in the 
TCNQ solids. Given the clear evidence for the ...1100... 
CO in the 1:2 TCNQ solids, the occurrence of the 
...1010... CO in TMTSF therefore appears questionable. 
Specifically, previous parameterizations of (TMTSF) 2 X 
have all put U/to < 6, which makes V c /to ~ 2.8. Taken 
together with the condition V < U/2, this makes the 
...1010.... order very unlikely. 

Direct evidence for the ...1100... order in (TMTSF) 2 X 
and (TMTTF) 2 Br comes from the determination that 
(a) there exists a coexisting CDW-SDW in these sys- 
tems below TsdWi an d (blithe periodicities of these 
density waves are both 2k F B0. Since the ...1010... CO 
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would give a 4kp charge periodicity, this can then im- 
mediately be ruled out. The coexisting CDW-SDW has 
been demonstrated in our previous_quasi-2D calculation 
within the Hamiltonian of Eq. (]2j)Ea, where we demon- 
strated that for small but nonzero interchain coupling 
the proper description of the ground state is a bond- 
charge-spin density wave (BCSDW), with the CO being 
...1100... 



C. (TMTTF) 2 X 

In a series of recent papers, Nad et al. have found 
evidence for CO in various (T ^ ^jffiT-E)*X materials, from 
dielectric permittivity studiesc3l3liJii3. Importantly, the 
CO is found both with centrosymmetric anions X = PFg, 
AsFg, SbFg etc. and with the asympetric anion Re04. 
Observation of ferroelectric behaviorlla suggests removal 
of all symmetry elements, and based upon this it has been 
suggested that the CO pattern here is of the ...1010... 
type&jEj and is accompanied by anion motion. Indepen- 
dent and direct verifications <ACO in X = PF 6 and AsF 6 
have come from NMR studiesO. However, while the pres- 
ence of some form of CO in (TMTTF) 2 X is undisputed, 
we believe there is not yet enough evidence to say defini- 
tively whether the pattern of CO is ...1010... or ...1100... 
Below we review the evidence from several different per- 
spectives and suggest a possible resolution. 

(i) The very high temperatures at which the CO is ob- 
served in some of these CTS (T co = 154 K and 227.5 K 
in X = SbF6 and Re04, respectively) indeed sugfiests-the 
...1010... order, since long range ...1100... orderOi§Ej is 
associated with the SP transition, which usually occurs at 
much lower temperatures. One possibility (not yet thor- 
oughly explored, let alone established) that would allow 
higher temperature CO within the ...1100... scenario is 
that in at least some of the (TMTTF) 2 X the CO is asso- 
ciated with pre-transitional structural fluctuations from 
the SP transition. Theoretically, it has been claimed that 
such fluctuations of the SP order may become visible at 
temperatures as high as 3-^- 4 T$p, where Tsp is the 
SP transition temperatureLa. Experimentally, structural 
fluctuations associated with the SP transitions have been 
observed in the related SP materials (BCPTTF) 2 PF 6 
and (BCPTTF) 2 AsF 6 at temperatures as high as 100 K 
and 120 K, respectively, even though their SP transition 
temperatures are T$p = 37K and 34 K, respectivelyE3E3. 
It is then intriguing that in (TMTTF) 2 PF 6 scattering at 
2kp in X-ray experiments begins to appear at neatly the 
same temperaturdfl at which CO appears in NMREJ. 

(ii) NMR experiments lead to estimates of a lower limit 
of 0.5 for AnEll, which also supports the ...1010... CO, 
since within our finite size calculations we have never 
found such a large An within the ID ...1100... phase. 
While the large An would be in agreement with the 
...1010... order, the difficult of determining An unam- 
biguously from these experiments is shown by the re- 
sults similar NMR experiments on the inorganic mate- 



rial NaV 2 05, in which also CO involving the vanadium 
sites is observed and in which the pattern of the CO 
remains controversial after considerable study. These ex- 
periments have not led to An values that are consistent: 
widely different An, ranging anywhere from 0.1 iitQ [ -Q.8 
have been reported in the literature in this caseE30Lj. 

(iii) Recent NMR experiments on X=AsF@ under pres- 
sure also find that the CO is very sensitive to pres- 
surej-apd is suppressed for pressures greater than ~0.15 
GPaEl The great sensitivity to pressure raises further 
questions regarding the CO in (TMTTF) 2 X: if the CO 
in (TMTTF) 2 X is of the ...1010... pattern and hence 
driven by large V, one would expect that increased tx/to 
from the application of pressure would have little effect 
on the CO. Importantly, studies of two- leg ^-filled lad- 
ders find V C (U) in the ladder to be quite similar to the 
value in ID for small t±/t(p3. 

(iv) It has been suggested that the MI transition at 
higher temperatures is not driven by intrastack Akp 
bond dimerization (which is unaccompanied by CO), but 
by the anion potentially. Ascribing the metal-insulator 
transition to the anion potential, however, is also prob- 
lematic. For instance, if the anion potential leads to a 
4kp- bond dimerization, there is further increase-in V C (U), 
and the ...1010... CO becomes even less likelyo. On the 
other hand, if the anion potential led to the ...1010... CO 
directly, evidence for the CO should have been found al- 
ready at the metal-insulator transition. As regards the 
widely differing Tcoj it is possible that different anion 
sizes lead to different interchain couplings, which in turn 
affect the V C (U). 

(v) Finally, why is Tco is so low in some of these sys- 
tems (Tcofi 65 K and 100 K in X = PF 6 and AsF 6 , 
respectively^)? The charge localization/metal-insulator 
transition T M i in (TMTTF) 2 PF 6 occurs at > 200 K 
above the CO temperature, while (as explained in Sec- 



tion 



IV| ) for the ...1010... model one expects Tco = Tj 



MI- 



A possible resolution to (i) — (v) is that the higher tem- 
perature CO found in (TMTTF) 2 X is indeed ...1010..., 
while the lower temperature SP state has ...1100... CO. 
Transitions at finite temperature are determined not by 
the ground state energy but by the free energy. At tem- 
peratures T < Tmi, the excitations of a strongly corre- 
lated system are predominantly spin excitations. Because 
of the greater multiplicities of high spin states, the free 
energy at high T (but below Tmi) is then dominated 
by high spin states. Consider now the fully ferromag- 
netic state, which is clearly described within the spinless 
fermion Hamiltonian limit of Eq. (Q). Thus V c in this 
state remains 2to even for finite U and is smaller than 
the V c for the ground state with total spin S = 0. For 
similar reasons we expect for finite U, V c to increase pro- 
gressively as the total spin S decreases. Therefore, for 
a system in which V is slightly lower than V c (S = 0) 
but higher than V c (S = S max ) we see that at high tem- 
peratures, where high spin excitations are thermally ac- 
cessible, the CO pattern can be ...1010... At lower tem- 
peratures where the high spin states become thermally 
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inaccessible and the free energy is dominated by low spin 
states there can be a switching to the ...1100... CO. We 
are current exploring this scenario in detail. 

Finally, we suggest two experiments that can settle 
the issue of CO in (TMTTF) 2 X. First, X-ray or neutron 
diffraction experiments can probe the actual bond distor- 
tions below the SP transition: with sufficient sensitivity, 
they should be able to distinguish the SW'SW distor- 
tions associated with the ...f fOO... CO from the SSWW 
distortions associated with the ...1010... CO. Second, as 



shown in Section HI, we predict three different charges 
below Tsp for the case of the ...1010... CO (rather than 
the two charges associated with the ...1100... state). If 
NMR experiments can be extended to this temperature 
region, further splitting of the NMR lines is predicted. 
Unfortunately, our results also suggest that the charge 
difference between the two '0's is rather small, and an 
extremely sensitive probe might be necessary to deter- 
mine the new splitting. 



VI. CONCLUSION 

Our conclusions can be summarized as follows. First, 
formation of the Wigner crystal-like ...1010... CO re- 
quires V within the ID extended Hubbard model that is 
much larger than that predicted within mean field the- 
ory. Further, mean field theory predicts the incorrect 
behavior of V C (U) as a function of U . Accurate many- 
body calculations show that for finite U, V C (U) > 2\to\ 
and increases with decreasing U. 

Second, we have given a complete theory of the SP 
transition that can occur in the correlated ID ^-filled 
band. The SP phase with the ...1100... CO is a BCDW, 
with two different charges and three different bonds, the 
bond distortion pattern being SW'SW. In contrast, the 
SP phase in the ...1010... CO has two kinds of bonds 
and three distinct charges; the bond distortion pattern 
here is SSWW, and this distortion makes the charges 
on the sites labeled '0' unequal. Several distinct types 
of experiments have confirmed the existence of two kinds 
of charges and three kinds of bonds in a number of 1:2 
anionic CTS. In (TMTSF) 2 X, experiments find coexist- 
ing CDW and SDW with the same periodicity, which as 
discussed above precludes the ...1010... CO. The situ- 
ation in the (TMTTF) 2 X is less clear, as experiments 
showing the existence of CO do not directly measure the 
specific CO pattern. We know of no experiments to date 
that have confirmed the existence of the three types of 



charges and two types of bonds associated with the SP 
phase of the ...1010... We believe that the NMR and X- 
ray experiments that could confirm this should have high 
priority. 

Third, within the strictly ID models, the two charge 
orders predict different scenarios for the temperature de- 
pendence of the CO. For the ...1100... CO, T MI > T co = 
T SP , whereas for the ...1010... CO, T MI = T co > T SP . 
Neither of these results is in agreement with all the exper- 
iments, which in general show Tmi > Tco > T$p. We 
believe that weakly 2D inter-chain couplings, which ex- 
plain the transition from CDWfSP to CDW-SDW ground 
states in some of the materialsEa, may also be the key to 
understanding the ordering of Tmi, Tco, an d T$p in the 
different CTS. Support for this belief comes from recent 
results on the strongly 2D 0-(BEDT-TTF) 2 X materials, 
in which Tmi occurs simultaneously with T ap ., fo llowed 
by a lower temperature spin-gap transitionHo. This 
sequence of transition temperatures, as well as the 2D 
CO and spin-gap formation, can be understood within 
a 2D, six-fold coordinated lattice model with dominant 
...1100... CO ordeiM 

Fourth, our calculations suggest several problems for 
further study. Apart from the experiments seeking to 
confirm a "three-charge, two bond" 4kp -CDW-SP state, 
these include theoretical studies providing a full expla- 
nation of the sequence of MI, CO, and SP transitions 
as the temperature is decreased, an accurate determina- 
tion of the phase diagram in the e-ph coupling constant 
plane in the thermodynamic limit, and a quantitative un- 
derstanding of the role of broken symmetries associated 
with higher spin states at finite temperature. We are 
currently investigating these theoretical issues. 

Finally, during preparation of this paper we became 
aware of the recent work of Shibata et al.E3, who investi- 
gated the purely electronic extended Hubbard model in 
ID. Although they did not discuss the ...1100... charge 
ordering, their results for V c are completely consistent 
with ours in Section EJ. They have also shown that dimer- 
ization increases V c . 
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